Dynamical chiral symmetry breaking in QED3 at finite density and impurity potential 
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We study the effects of finite chemical potential and impurity scattering on dynamical fermion 
mass generation in (2+l)-dimensional quantum electrodynamics. In any realistic systems, these 
effects usually can not be neglected. The longitudinal component of gauge field develops a finite 
static length produced by chemical potential and impurity scattering, while the transverse com- 
ponent remains long-ranged because of the gauge invariance. Another important consequence of 
impurity scattering is that the fermions have a finite damping rate, which reduces their lifetime 
staying in a definite quantum state. By solving the Dyson-Schwinger equation for fermion mass 
function, it is found that these effects lead to strong suppression of the critical fermion flavor 
and the dynamical fermion mass in the symmetry broken phase. 
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I. INTRODUCTION 

The dynamical chiral symmetry breaking (DCSB) 
in (2-|-l)-dimensional quantum electrodynamics (QED3) 
has been investigated intensively for more than twenty 
years [ll-IT^. On the one hand, these investigations may 
help to gain deeper understanding of DCSB in QCD. On 
the other hand, this non-perturbative phenomenon gives 
an elegant field theoretic description for the formation 
of long-range antiferromagnetic order in two-dimensional 
quantum Heisenberg antifcrromagnet 

Using the Dyson-Schwinger (DS) equation approach, 
Appelquist et al. first fomid that DCSB caused by 
fermion mass generation takes place only when the 
fermion flavor N is below some critical value Nc Q. 
When N > Nc, the fermions remain massless and the 
chiral symmetry is preserved. Although there was some 
controversy about the existence and precise value of Nc 
in the literature, most analytical and numerical compu- 
tations confirmed that TV^ « 3.5 0-11, i-IHl H • 

The DCSB is realized by forming fermion-anti-fermion 
pairs mediated by strong gauge interaction. In order 
to trigger this low-energy phenomenon, an essential re- 
quirement is that the gauge interaction between fermions 
should be sufficiently long-ranged. If the gauge field de- 
velops an effective screening length, then the gauge in- 
teraction is significantly weakened. We have considered 
the effect of a finite gauge boson mass rria on DCSB pre- 
viously [l7| , and found that the increasing rua leads to a 
significant suppression of DCSB. 

When the action contains only massless Dirac fermions 
and non-compact gauge field, the gauge interaction is al- 
ways long-ranged since the gauge invariance ensures that 
no photon mass term can be generated by the radial cor- 
rection from fermions. If the gauge field also couples to 
an additional Abelian Higgs model, then the photon can 
acquire a finite mass via the Anderson-Higgs mechanism 
in the superconducting state [13] . In the context of high 
temperature superconductor, the suppression of DCSB 
by photon mass ttIq can be used to qualitatively under- 
stand the competition between antiferromagnetism and 



superconductivity, which is one of the most prominent 
phenomena in high temperature superconductors [l7| . 

Apart from Anderson-Higgs mechanism, the gauge 
field may also be screened by other physical effects. For 
example, when the fermion density is nonzero, the chemi- 
cal potential /i should be explicitly taken into account. In 
addition, the fermions may be scattered by some impu- 
rity atom (defect or imperfection) in realistic interacting 
system described by QED3 . Although generally chemical 
potential and impurity scattering have distinct influences 
on the physical properties of the interacting fermion sys- 
tem, there is a common feature that they can both induce 
a finite density of states at the Fermi level. This screens 
the time component of the gauge field, which then be- 
comes short-ranged. However, the transverse gauge field 
remains long-ranged in the sense that the static screen- 
ing length still diverges as at zero chemical potential in 
the clean limit (" clean" means there is completely no im- 
purity scattering). In principle, the long-range nature of 
transverse gauge field is guaranteed by the gauge invari- 
ance. However, though the transverse part is unscreened, 
the effective gauge interaction strength is reduced by the 
fermion density caused by chemical potential and/or im- 
purity scattering process, which certainly affects the crit- 
ical flavor Nc as well as the dynamical fermion mass. Be- 
sides generating flnite fermion density, the impurity scat- 
tering leads to the damping of fermionic quantum states 
by producing a flnite scattering rate Fq. Such damping 
effect reduces the lifetime of a fermion staying at a defl- 
nite quantum state, specifled by such quantum numbers 
as momentum and/or energy. This will also have impor- 
tant effects on DCSB. 

While there appeared several papers discussing the ef- 
fect of chemical potential on DCSB in QED3 [l^, to 
our knowledge the effect of impurity scattering has never 
been considered in previous work. Indeed, in any real- 
istic applications of QED3 to condensed matter physics, 
such as high temperature superconductor, the chemical 
potential and impurity scattering are usually important 
and thus should not been ignored. The purpose of this 
paper is to study the dependence of DCSB on flnite den- 
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sity and impurity scattering. 

The paper is organized as follows. In Sec. [Ill we set up 
the Lagrangian and then discuss the screening effect of 
gauge interaction by including chemical potential and im- 
purity scattering rate into the vacuum polarization func- 
tion. We then solve the Dyson-Schwinger equation and 
present the critical flavor and fermion mass function in 
Sec. mil The discussion is presented in Sec. |lVl 



II. 



DYSON-SCHWINGER EQUATION IN THE 
PRESENCE OF ^ AND To 



The Lagrangian for (2+l)-dimensional QED with N 
flavors of massless fermions is given by 



N 



C = -lF,,F'^-' + yyjt^^-e4^)^a. (1) 



=1 



In (2+l)-dimensional space-time, the lowest rank spino- 
rial representation is two-component spinor whose 2x2 
representation may be chosen as the Pauli matrices 
7^ = ((72, iCTs, i(Ti). However, it is impossible to define 
a 2 X 2 matrix that anticommutes with all these ma- 
trices. Therefore, there is no chiral symmetry in this 
representation. The fermion can be described by a four- 
component spinor field Tp, whose conjugate spinor field 
being defined as ip — tp'^jo The 4x4 7- matrices 
can be defined as 7^ = (cts, ici, j(T2) (8) c^, satisfying the 
standard Clifford algebra {7^^,71/} = 217^1, with metric 
5^,^ = diag(l, —1, —1). It is easy to verify that there are 
two 4x4 matrices 



73 = « 



which anticommute with all 7^. The massless Lagrangian 
(Eq. ([1])) preserves a continuous U(2N) chiral symmetry 















75 = M 





'"ijj. The mass term generated by fermion- 



anti-fermion pairing will break this global chiral sym- 
metry dynamically to subgroup U{N) x U{N). In the 
following we consider a general large N and perform the 
1 /N expansion. For convenience, we work in units where 
h = ks = ^■ 

In the Euclidian space, the full gauge field propagator 
Diivisi) is given by the equation 

D-'Aq)^Dfy\q)+Ii,,{q)^ (2) 

with the free photon propagator being 

(3) 



in the Landau gauge. To the leading order of expan- 
sion, the one-loop contribution to vacuum polarization 
tensor 11^^ is 



d^k Tr[7,j7,(^ + ^)] 



(4) 



where a = Ne^ . The QED3 can be treated using the \/N 
expansion, with the product a = Ne^ being fixed when 
N ^ 00 and 0. This tensor can also be written 

as n^^(g) = (g^^ - ^||^)n(g2) according to the gauge 
invariance. It is easy to find the polarization function: 
nl?) = Now the propagator of gauge field has the 
form 



(5) 



In order to study the dynamical fermion mass genera- 
tion, one can write the following DS equation for fermions 
propagator 



S-/{p)^sf-\p)^e^ 



(Pk 

-^—^-1^SF{k)D^,^{p-k)T^{k,p). 



(6) 

To the leading order in 1/iV expansion, the vertex func- 
tion is replaced by the bare matrix 7^. The inverse 
full propagator of fermion is 



Sf{p)-^ ^ii>+Y.{p), 



(7) 



where the wave-function rcnormalization is neglected. 
Taking trace on both sides of the DS equation, we get 
an integral equation 



(fk 



-Tr[7^i?^,(p-fc)7,]. (8) 



4 J (27r)3fc2-hS2(fc) 
Using the gauge field propagator Eq. (O, we have 
(fk S(fc) 1 



2a 



N J {2nfk^ + j:^{k)(p-k)^ + 2(P_l^' 



(9) 



which was first obtained by Appelquist et al. Q. They 
found that DCSB takes place only when the fermion fla- 
vor is less than a critical value Nc = 32/7r^. Motivated by 
this interesting prediction, a great many attentions have 
been paid on this issue. After taking into account the 
high order corrections to wave function rcnormalization, 
the critical flavor was found to change to Nc = 128/37r^ 
Q . Pennington et al [3] used a more careful truncation of 
the fermion DS equation and found that DCSB can occur 
for all values of N . When the DS equations of fermion 
self-energy is coupled self-consistcntly to those of gauge 
field propagator. Maris Q showed that iVc ~ 3.3. In a 
recent publication, Fischer et al. [ll| found that Nc — ^ 
after detailed analysis of vertex corrections. Some nu- 
merical simulations on lattice QED3 claim that there is 
no DCSB for iV > 2 12]. However, Gusynin et al. found 
that the absence of DCSB can be attributed to the large 
infrared cutoff used in lattice studies and the smallness of 
the generated mass scale [ll] . In summary, most analyt- 
ical and numerical computations seem to agree that the 
critical fermion flavor should be Nc ~ 3.5 [sl-lsl. l9l-[TTi [l^ . 
close to the original value obtained by Appelquist et al. 
within the lowest order of 1 /N expansion. 
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The above investigation and result are valid at zero 
chemical potential in clean fermion system. Till now, lit- 
tle attention has been paid to the case with finite fermion 
density and finite impurity potential. Since the chemi- 
cal potential and impurity scattering are generally very 
important in realistic applications, it is important to ex- 
amine their effects on DCSB. The aim of this paper is to 
study this problem. 

First of all, at nonzero chemical potential the fermions 
have a finite density at the Fermi level, which screens 
the temporal component of gauge field and weakens the 
gauge interaction. The influence of chemical potential 
on DCSB was previously discussed by Feng et al. [l8| . 
They derived an explicit equation to include the chemical 
potential into the DS equation. One crucial assumption 
in their work is the neglecting of chemical potential in the 
vacuum polarization function. The numerical results of 
Feng et al. show that the chemical potential leads to only 
insignificant change of N,. and fermion mass. We think 
that the screening effect caused by chemical potential is 
very important and hence will pay special attention to the 
screening of gauge field induced by chemical potential by 
studying the vacuum polarization. As will be shown in 
the context, the chemical potential suppresses both Nc 
and fermion mass strongly. 

The effect of impurity scattering is more complicated 
than chemical potential. Generally, the scattering of 
fermions by impurity potential has two important effects. 
First, it generates a finite density of states at low energy, 
which also screens the gauge field. This effect is expected 
to weaken the gauge interaction, analogous to the role 
played by the chemical potential. Second, the impurity 
scattering produces a finite damping of fermion quan- 
tum states and thus reduces the time for massless Dirac 
fermions to interact with their anti-particles. These two 
effects are both very important. 

The behavior of massless Dirac fermions in a random 
potential is of great interest in the context of condensed 
matter physics since the low-energy properties of some 
many-body systems, such as rf-wave high temperature 
superconductor and graphene, are largely controlled by 
the interaction of Dirac fermion with impurities. Unfor- 
tunately, at present the problem of random Dirac fermion 
has not been fully understood even when there is no di- 
rect interaction between fermions. When the impurity 
scattering and gauge interaction are both important, the 
problem is basically out of theoretical control. If we con- 
sider a single impurity atom, then the impurity scattering 
can be treated by the self-consistent Born approximation. 
Within this approximation, the retarded fermion self- 
energy function develops a finite imaginary part, which is 
usually represented by a constant scattering rate Fq. To 
study the problem about impurity scattering, it is most 
convenient to work in the Matsubara formalism and write 
the fermion propagator as 



where the frequency is iuj„ 



;(2n+l)£ 

13 



with /3 = y . Once 



S'_F(iw„,p) 



1 



the impurity scattering rate is taken into account, the 
fermion frequency should be replaced by 



iw„ ^ iw„ + iFoSgn(w„). 



(11) 



The fermion damping effect can be intuitively understood 
as follows. After analytical continuation. 



iuJn — > o; + «(5, 
the retarded fermion propagator has the form 

1 



(w + iFo)7o - 7 ■ P' 
After Fourier transformation, it becomes 



(12) 



(13) 



(14) 



in real space. Apparently, the parameter Fq measures 
the decaying rate of the fermionic state characterized by 
such quantum numbers as (w, p), which is known as the 
Landau damping effect. Starting from the propagator 
Eq. (1131) . various physical quantities can be calculated 
and compared with experiments. Intuitively, the damp- 
ing effect is at variance with DCSB since a fermion may 
be scattered into another state before it combines with 
its anti-particle to form a stable pair. 

Before we set up the DS equation for fermion self- 
energy, we need first to calculate the photon propaga- 
tor and discuss the screening effect caused by chemical 
potential and impurity scattering within the Matsubara 
formalism. With finite chemical potential /i and finite 
damping rate Fq, the fermion propagator reads 



1 



{iLOn + iroSgn(iLj„) - ^)7o - 7 • P 



(15) 



lUJnlO - 7 • P 



(10) 



The energy shift caused by chemical potential and the 
damping effect caused by impurity scattering are both 
reflected in this propagator. The screening of gauge in- 
teraction induced by them can be seen from the corre- 
sponding vacuum polarization functions. We will use this 
propagator to calculate the vacuum polarization func- 
tions and to construct the DS equation for fermion mass. 

The inverse photon propagator for frequency = 2^ 
and spatial momentum |q| is given by 

A;J(go,q,/3) = AjS"\go,q,/3) + n^.(go,q,/3), (16) 

where the IS.^^} is the free photon propagator. Here, we 
use different symbols and A^,y to denote the gauge 
boson propagator at zero temperature and finite temper- 
ature, respectively. Taking advantage of the transverse 
condition, qfjl^j,y{q) — 0, the vacuum polarization tensor 
defined by Eq. ^ can be decomposed in terms of two 
independent transverse tensors 'lO], 

n^i^((7o,q,/3) = ^A{qQ,(i,l3)Af,^ -i-nB(go,q,/3)s^i/,(i7) 
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where 



larization tensor II^i^ by the following expressions 



(19) 



They are orthogonal and related by the relationship 

A 4- R — A — 



The functions n^(q'o,q, /3) and 113(50, q, /?) are related 
to the temporal and spatial components of vacuum po- 



2 2 
= ^Hoo, Hb = Tlu - ^Hoo. (21) 



Now the full finite temperature photon propagator 
A^,y(go, q, /?) can be written as 



(20) A^^(9o,q,/3) 



A 



B 



g2 + n^(go,q,/3) g2 + nB(go,q,/3)' 

(22) 

At finite temperature, the fermion contribution to vac- 
uum polarization functions should be 



J 



n— — 00 



■7 1 - s 



n— — 00 



d^L 2l^l,y + (1 - 2x){l^q^ + q^_,l^ - q ■ IS^^) + 2x{l - x){q^5^^ - q^q^) 



(2^)2 



[P + x{l - x)q^ 



)^P + x{l-x)q^' 



(23) 



where 9 = (90, q), Q = |q|, 90 = for gauge boson and 
k ~ (fco,k), K — |k|, fco = fgj-jjijon. Here, a 

new momentum variable is defined hy I — k + xq with 

l^{lo,l), L=\l\, = ^(n + xm + i). (24) 

In the presence of impurity scattering rate Fq and chem- 
ical potential /i, the variable Iq should be replaced by 

2tt 1 /3 /3 

'o(") = — (ri + a;m + - + i—^ + — roSgn(w„)). (25) 

Now the spatial and temporal component of polarization 
tensor n^i, can be expressed as 



n,(go,q,/3) = y^ dxj ^ 
x[2x(l - x){q^S,j - g,(7,)52 - (1 - 2x)5,,go^*](26) 



noo('7o,q, fi) 



4a 



dx 



d^L 



P Jo "-^J (27r)2 
X [Si - 2(l2 -I- x{l - x)ql)S2 + (1 - 2a;)go5*] (27) 
where, following Ref. |19j, we defined 

1 



n— — c 

00 

= E 



[;2(n)+i2+a;(l_2;)g2 

^o("-) 

[ll{n) + + x{\ ^ x)q' 



(28) 
2- (29) 



It is not easy to calculate Eq. (|26p and Eq. ([77]) ana- 
lytically. Nevertheless, within the widely used instanta- 
neous approximation q^ — Q, the integration can be per- 
formed by the methods presented in Ref. fio'l and 
After tedious but straightforward computation, we ob- 
tain the following expressions for polarization functions: 



nA(q,/3) 



X In 



nB(q,/3) 



2a 

2 ( cosh 

Aa 

T 



dx 



10'"2^r + 1 



2'KXq{x) 



10^''^ Xr + 1 
1 

dx Xq{x) 





cosh 



2nX„, 



10'" 2;^^, _|. I 



sinh 






10i"2Xr + l 


cosh 


2TrXq{x) 


+ cosh 


27rX^ 


10i»2Xr + l 


10'"2Xr + l 



(30) 



where Xq{x) = £^xil-x)q, Xr = J^Eo, X 



In either of the limits /3 ^ 00 (T 0) or q — > 00, both 
the above functions reduce to aq/8 and then the photon 
propagator reduces to the zero-temperature result (Eq. 
([5])). However, in the limit q — 0, Hyi(q, /3) becomes a 
function of temperature, chemical potential, and impu- 
rity scattering rate, Moo{f3 , Xr , X fj,) . In the zero energy 
and zero momentum limit, the photon propagator is 



A^.(go =0,q^0,/3) 



A 



B 



q2 + Mi, 



(31) 



q2 + 

Apparently, the temporal component of the propagator 



5 



now acquires an effective mass and becomes 



Moo = 



1 2a 



(10i"2j^r + l)ln 



2 cosh 



10'" ^Xr + 1' 



(32) 

The mass appearing in the longitudinal photon implies 
that the electric field acquires a static screening length 
determined by the chemical potential and impurity scat- 
tering rate. The transverse photon, however, remains 
massless and hence the corresponding magnetic field is 
still long-ranged, albeit dynamically screened. 

We now study the problem of DCSB at finite temper- 
ature. As in the case of zero temperature, this problem 
can be studied by the DS equation approach. The DS 
equation for fermion propagator Sp at finite tempera- 
ture is given by 



E 



2'^A"S'F(^o,k,/3)r^A^i,(9o,q,/3), 

(33) 



where A{po,P, is the wave- function renormalization. 
Substituting Eq. ^ and Eq. JM]) into Eq. JM]) and 
then taking trace on both sides of this equation, we get 
the following couple of closed integral equations 



^(po,P,/3) = 1 



-y 

4/3 ^ 



(27r) 



rTr 



S(po,P,/3) = -^^ 



(2^ 



Tr [7^S'Fr^A^i,] 



(35) 



Here, A^j^ is the photon propagator as defined by Eq. 
(IT51) and the is the full vertex function. To treat this 
equation, a number of approximations should be made. 

At present, there is no well-controlled way to choose 
the vertex function T^. In order to satisfy the Ward- 
Takahashi identity. Maris et al. Q studied several dif- 
ferent Ans'dtze for the full vertex function and compared 
the results. For the sake of simplicity, they assumed that 
= f{A{p),A{k),A{p ~ k))-f^, with A{p) = 1 corre- 
sponding to the bare vertex. They found that, within a 
range of vertex functions, the critical behavior of DCSB is 
almost independent of the precise form of /-function and 
the bare vertex is actually a good approximation. Once 
Ti/ = 7i/ is assumed, then the Ward-Takahashi identity 
requires that the wave function renormalization should 
be A{p) = 1. 

In the present case, it is much more difficult to solve 
the DS equations Eq. (I55|) with finite chemical potential 
and impurity scattering than in the clean system. Unlike 



the zero temperature case, the particle energy and mo- 
mentum should be treated separately in the finite tem- 
perature perturbation theory. Now, the energy becomes 
discrete frequencies proportional to the temperature T. 
To simplify the theoretical and numerical analysis, we 
keep only the leading order of 1 /N expansion by replac- 
ing the vertex function T^, by 7^ and taking A{p) = 1. 
The DS equation of fermion Eq. p4p can now be simpli- 
fied as 



Sp^{po,P,f3) = [(iw„ + iroSgn(cj„) - ^)7o - 7 • p] 



+S(po,P,/3). 



(36) 



The fermion mass function E(po,P,/3) satisfies the fol- 
lowing integral equation 

nPo,p,f3)^jY. J (^Tr[7^7,A^,(go,q,/3)] 

^ S(fco,k,/3) 

^ 4 [{un + Tosgniun) + j/i)^ + k2 + S2(fco, k, /?)] 



a f3 
iV4^ 



E 



d^k 



where Tu is the fermion-photon vertex and q = p ^ k. 
The full fermion propagator can be written as 

5pi(po,P,/3) = ^(po,P,/3)4°^"' + 5^(Po,P,/3), (34) where 



(2^)2 
E(fco,k,/3) 



A^^(go,q,/3) 



[{{n + i) + Xrsgniuj„) + iX^)^ + M^] 



27r' 



(37) 



Xt^-^, M = XT\/k2-f S2(k,/3). (38) 



Note that chemical potential and impurity scattering rate 
both appear in two places: in the occupation number 
and in the polarization function. The energy shift in- 
duced by chemical potential and the fermion damping ef- 
fect induced by impurity scattering are both represented 
in the former place. The screening of gauge interaction 
caused by chemical potential and impurity scattering is 
reflected in the vacuum polarization function. In order 
to perform the frequency summation in Eq. (j37p . we uti- 
lize the approximation I](po,P, /3) — S(po = 0,p,/3). At 
t^rn = 90 = 0, the summation over a;„ in Eq. p7|) yields 
the equation 



S(p,/3) 



13a 
4iV7r2 



rf2k I](k,/3) 

W) 



M 



App(go = 0,q,/3) 



i;[-+Xr±iX^ + iM] 



(39) 



Substituting the full photon propagator Eq. into the 
mass equation and use the same approximative method 
as that presented in Ref. p^l, we have 



S(p,/3) 
Pa f d'^k 



4iV7r2 J {2tt) 
S(k,/3) ^ 
M 



(72+n^(q,/3) g2 + nB(q,^) 



i^[-+Xr±iX^ + iM] 



6 



f3a 



ANtt^ J (27r)2 
•E(k,/3)^ 



1 



(72+n^(q,/3) q2 



M 



tanh 



10i"2Xr + 1 



ni3(q,/3). 

(40) 



If this nonlinear integral equation develops a nontrivial 
solution, then the massless fermion acquires a finite mass. 



III. SOLUTION OF THE DYSON-SCHWINGER 
EQUATION 

The nonlinear integral equation can be solved by the 
bifurcation theory and parameter imbedding method 
[21I [2^ . The basic idea and detailed computation pro- 
cedures are presented in previous papers [13, HH US] • To 
determine the bifurcation point that separates the chiral 
symmetric phase and symmetry broken phase, we need to 
find the eigenvalues of the associated linearized equation. 
Taking the Frechet derivative of the nonlinear integral 
equation Eq. pOj) . we have 



S(p,/3) 

l3a 



1 



1 



47V7r2 ./ (27r)2 

d ri](k,/3)7r 



g2 + n^(q,/3) g2 + nB(q,/3) 



xS(k,/3)^^ 

(3 r rf2it 

ANtt^ J (27r)2 
. S(k,/3)^ 



M 



2 
1 



tanh 



Tr{M ± X^, 



10 



In 2 



Xr + 1 
1 



Xrlkl 2 



tanh 



7r(XT\k\±X^) 
10^'^^ Xr + 1 



ni3(q. 



(41) 



To facilitate numerical calculation, we divide the mo- 
menta, temperature, fermion mass, chemical potential 
and impurity scattering rate by parameter a and make 
them dimensionless. The UV cutoff for momentum now is 
replaced by A/a, which should be properly chosen to en- 
sure the results insensitive to UV cutoff. As pointed out 
by Appelquist et al. [2, Q , the integral in the zero tem- 
perature DS equation damps rapidly for momenta greater 
than a, it is thus natural to set A/a ~ 1. This is also 
true for DS equation at finite temperature [7|. In the 
following we take the advantage of this fact and assume 
that A/a = 1. 

The dependence of critical flavor number Nc on chem- 
ical potential ^ and impurity scattering rate Tq are 
shown in Fig. [1] for a list of values of temperature 
T = 10-^,10-^^,10-^. It is easy to see that the criti- 
cal flavor Nc is a decreasing function of T, /i, and Fq, 
implying a strong suppression of DCSB by thermal fluc- 
tuation, fermion density, and impurity scattering. 

The fermion mass function S(p, /3) can be obtained af- 
ter solving the DS equation Eq. (|40|) using the straight- 
forward iteration method. The results are presented in 
Fig. [5] and Fig. [3] In Fig. ^ we show the re-scaled zero- 
momentum fermion mass S(0) for several different values 




FIG. 1: Dependence of Nc on n and To for different temper- 
atures T = 10-^ 10-'' and IQ-". 



of fermion flavor at a fixed temperature T — lO^^*', 
with /i = Fq = 0. Apparently, the fermion mass S(0) 
decreases rapidly as the fermion flavor A^ increases. This 
is easy to understand since the inverse of fermion flavor 
1/A^ serves as the effective coupling strength of gauge 
interaction. 

The dynamical fermion mass as a function of momen- 
tum is shown in Fig. [3] for A^ = 2 at different values of 
temperatures, chemical potential and impurity scattering 
rate. We notice that the dynamical fermion mass is sig- 
nificantly suppressed by the increasing chemical potential 
and impurity scattering rate. It also decreases with the 
increasing fermion momentum. 

We next would like to compare the results presented 
above with those in QED3 with finite photon mass gen- 



1.6 2.0 2.4 2.8 

N 



10"" 
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FIG. 2: Zero-momentum mass E(0) for different values of N 
at T = 10"^° and ^ = To = 0. 



erated by Anderson-Higgs mechanism. Different from 
the partial screening of gauge interaction, the Anderson- 
Higgs mechanism induces a complete static screening. In 
this case, the gauge boson acquires a physical mass by 
eating the massless Goldstone bosons, so that the prop- 
agator (Eq. ([22)) ) becomes 



A^^((7o,q,/3) = 



A 



g2 +nB(go,q,^) 



(42) 



If nia 7^ 0, both temporal and spatial components of 
gauge field develop a finite screening length, which is 
given by the gauge boson mass. Such screening length 
eliminates the contribution of small momenta to the DS 
equation Eq. (HT|) . It is expected that a large nia will 
prevent the DCSB. The dependence of critical flavor Nc 
on rua and Tq is shown in Fig. 4. The gauge boson mass 
and impurity scattering lead to similar suppression effect 
of DCSB. 



IV. SUMMARY AND DISCUSSION 

In summary, we studied the effects of finite chemi- 
cal potential and impurity scattering on dynamical mass 
generation in QED3. In the realistic applications of 
QED3 to condensed matter systems, these effects usu- 
ally can not be ignored. By solving the DS equation 
for fermion mass function, we found that both chemical 
potential and impurity scattering lead to strong reduc- 
tion of critical fermion flavor Nc and dynamical fermion 
mass. In reality, the chemical potential, impurity scat- 
tering, and even a finite gauge boson mass may coexist at 
the same time. When their effects are all important, the 
DCSB is completely suppressed. These results impose an 
important constraint on the applicability of QED3 in var- 
ious physical systems. For a system with large fermion 
density and impurity potential, it seems impossible that 
DCSB can indeed take place. 




FIG. 3: (a) Zero-momentum mass S(0) for different values 
of n and Fo at T = 10~*; (b) Dependence of mass on 
momentum p for different values of ^ and Fo at T = 10"''. 




FIG. 4: Dependence of Nc on ma and Fo at T = 10 



We now would like to briefly discuss the meaning of 
flnite chemical potential and the possible experimental 
study of its effect on DCSB. In condensed matter physics, 
the zero temperature ground state (vacuum) of a many- 
body system is characterized by the presence of a Fermi 
level which separates the fully occupied and fully empty 
states. In some condensed matter systems, such as high 
temperature superconductor and graphene, the valence 
band and conduction band touch at discrete Dirac points. 
When the Fermi level lies exactly at Dirac points, the 
chemical potential is zero and the low-energy excitations 
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are masslcss Dirac fcrmions. The Fermi level moves 
upwards (downwards) from the Dirac points once the 
fermion density increases (decreases). Now the chemi- 
cal potential is just the quantity that measures how the 
new Fermi level is far from the original Dirac points. At 
zero chemical potential, the pairs are formed by fermions 
slightly above Dirac points and anti-fermions (holes in 
condensed matter terminology) slightly below the Dirac 
points. Once the fermions become massive, they are con- 
fined. At finite chemical potential, the pairs are formed 
by fermions slightly above the finite Fermi surface and 
anti-fermions slightly below Fermi surface. The massive 
fermions are also confined by the gauge force. 

In realistic condensed matter systems, the fermion den- 
sity or chemical potential can be continually turned by 
chemical doping or gate voltage. We first have a number 
of samples, each with different chemical potential, and 
let them stay at some finite temperature where DCSB is 
completely prevented by thermal fluctuations. We then 
cool down these samples and study whether DCSB in- 
deed take place at very low temperature by measuring 
some observable quantities, such as specific heat, suscep- 
tibility, and electric conductivity. Through this way, in 
principle the effect of finite chemical potential on the fate 



of DCSB can be experimentally studied. 

We emphasize that our treatment on the random po- 
tential applies only to the case of weak impurity poten- 
tial where the self-consistent Born approximation is valid. 
When the random potential caused by impurities takes 
the form of Gaussian white noise, such treatment might 
be questionable. For this kind of impurities, the random 
potential has to be averaged by performing proper func- 
tional integration, which yields an effective four-fermion 
interaction in the whole action. This new effective inter- 
action can contribute a new term to the DS equation. At 
present, it is unclear technically how to study the DCSB 
and the fermion damping effect in a unified framework. 
This problem is subject to further investigation. 
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